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The western spruce budworm is an important pest of 
several conifers in western North American forests. A 
model called the Budworm Dynamics Model has been 
constructed that represents the population dynamics 
and feeding. This model is linked to the Prognosis 
Model, an individual tree, distant-independent stand 
growth model widely used in western North America. 

In this paper, a foliage dynamics model is presented 
that translates the tree attributes carried in the Progno- 
sis Model into foliage biomass that serves as a sub- 
strate for budworms. A tree damage model is pre- 
sented that translates budworm-caused defoliation into p a p 
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INTRODUCTION 

The Western Spruce Budworm Modeling System is a linked set of 
models that represent forest stand growth and western spruce budworm 
(Choristoneura occidentalis Freeman) population dynamics. The interac- 
tions between the trees and budwonns are explicitly represented in the sys- 
tem. The Prognosis Model for Stand Development (Stage 1973; Wykoff and 
others 1982) is a major component of the system. Its tasks include predict- 
ing the growth, mortality, and establishment of trees in a stand under 
budworm-free conditions for one or more rotations of trees. A Budworm 
Dynamics Model (Crookston and others 1990; Sheehan and others 1989) 
simulates budworm population dynamics through a short period, usually 
10 years. The Budworm Dynamics Model represents all life stages of the 
insect, accounts for several causes of budworm mortality, simulates adult 
dispersal, and computes the amount of foliage consumed by larvae. 

The state variables that represent stands in the Prognosis Model are not 
the same as those used by the Budworm Dynamics Model. For example, 
the Prognosis Model does not normally compute the amount of foliage bio- 
mass on each sample tree it represents. This information is needed by the 
Budworm Dynamics Model. Furthermore, the Prognosis Model normally 
does not represent budworm-caused damage in its predictions of tree 
growth. This step is necessary to represent the interactions between bud- 
worm and forest dynamics. 

This paper describes model components that predict the amount of foliage 
available for budworms to eat and those that predict the impact budworm- 
caused damage has on tree growth, mortality, topkill, and regeneration 
success. It is a companion paper to Sheehan and others (1989) who de- 
scribed the structure and content of the Budworm Dynamics Model. A de- 
scription such as this one is an abstract of the actual model. The source 
computer code is the only complete statement of the model (see inside back 
cover for information on acquiring Budworm Model programs). 

A certain amount of background information is necessary to understand 
this paper. Consult Wykoff and others (1982) for information on the Prog- 
nosis Model. For background on the Budworm Modeling System, consult 
Sheehan and others (1989), Crookston and others (1990), and Crookston 
(1991). 

The design criteria that influenced the development of this model are re- 
viewed in the following section. The foliage dynamics model follows, along 
with available references to the ideas it contains and the scant data on 
which it is based. The behavior of the model is presented without comparing 
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it to actual data. Next are the models that translate defoliation into tree 
damage, followed by illustrations of these models' behaviors. 

Some of the model components are based on little hard evidence. Future 
research needs are presented that summarize important gaps in our 
knowledge. 

A comment about the measurement systems used in this paper is in or- 
der. The Prognosis Model reports output data in English units of acres, 
feet, inches, and so on. The Budworm Dynamics Model is a product of a 
joint effort between the United States and Canada and the units of meas- 
urement are all metric. The result is that both measurement systems are 
used in the combined model. In this paper, the units reported are those 
program users are likely to see in output reports. Conversions between 
measurement systems are provided. 

DESIGN CRITERIA 

The Budworm Model was designed to forecast budwonn dynamics, tree 
damage, and subsequent stand growth and development. Forest cultural 
activities can be done that mitigate losses due to budworm. The models 
are used to simulate the results of management actions, and the simula- 
tion results are used to evaluate the benefits of the management actions. 
Therefore, one design criterion for the foliage models is that they may be 
made sensitive to the consequences of forest management practices. 

Table 1 contains a list of species represented by the Budworm Modeling 
System. Host species represented by the model are indicated in the table. 
The nonhost species are listed because estimates of total foliage biomass of 
nonhost also are required to represent some aspects of budworm dynamics. 

Budworm feeding and growth rates depend partly on the age-class of foli- 
age being consumed. The amount of new foliage, 1-year-old, 2-year-old, 
and remaining foliage needs to be known to the Budworm Model. Three 
tree-size classes for each host species are used to represent a forest stand 
(table 2). The crown of each tree class is divided into thirds. 



Table 1 Species represented by the Budworm Modeling System 



Common name 


Host 
status 


Abbrevi- 
ation 


Scientific name 




Douglas-fir 
Engelmann spruce 


Host 
Host 


DF 
ES 


Pseudotsuga menziest'i (Mirb.) 
Pteea engelmannii Parry 


Franco 


grand fir 

lodgepole pine 
ponderosa pine 


Host 
Nonhost 
Nonhost 


GF 
LP 
PP 


Abies grandis (Dougl.) Lindl. 
Pinus contorta Dougi. 
Pinus ponderosa Laws. 




subalplne fir 
western hemlock 
western larch 


Host 
Nonhost 
Nonhost 


AF 
WH 
WL 


Abies lasiocarpa (Hook.) Nutt. 
Tsuga heterophylla (Raf.) Sarg 
Larix occidental's Nutt. 





western recfcedar 


Nonhost 


we 


Thuja plicataDonn. 




western white pine 
white fir 1 


Nonhost 
Host 


WP 
WF 


Pinus monticola Dougl. 
Abies conco/or(Gord. & Glend. 


) Lindl. 



m T l !h P reS0nted Sf the ! nland F mpire Variant of th9 Prognosis Model, The Budworm Model 
may be linked to other geographic variants of the Prognosis Model and they may represent this species. 



Table 2 Foliage cell (crown third) codes used by the Budworm 
Modeling System 



Tree height 
size class 1 


Tree class 
name 


Crown third 
name 


Foliage cell 
numeric code 


0-7 m 


small 


top 

middle 


1 
2 






bottom 


3 


7-1 4m 


medium 


top 

middle 


4 
5 






bottom 


6 


Over 14m 


large 


top 
middle 


7 
8 






bottom 


9 



7m =23 ft, 14m = 46 ft. 



Include the 
Prognosis Model 
for Stand 
Development 



The Prognosis Model and thus the state variables it contains, its logical 
flow, and program architecture was chosen as a component of the Budworm 
Modeling System regardless of any negative consequences that decision may 
have had on the development of other Budworm Model components. The over- 
whelming success, flexibility, and user base of the Prognosis Model made it 
an extremely attractive cornerstone on which to build. This context for the 
foliage and damage models influenced their development in the following 
ways: 

1, To predict growth, the Prognosis Model typically uses forward differ- 
ence equations, most of which were calibrated on 10-year-long increment 
data. The Budworm Model uses 1-year-long increments that contain even 
shorter increments when some budworm development stages are simu- 
lated. The structure of the foliage model needs to allow periodic (10-year) 
updating of tree and foliage profiles while allowing annual updating of the 
foliage as it is influenced by defoliation. The structure of the damage mod- 
els must allow for consistent periodic updating of the Prognosis Model 
sample trees, even though the components of the damage model were cali- 
brated for specific periods that differ from the Prognosis Model and from 
each other. 

2. The Budworm Model operates on three strictly defined tree-size 
classes, and the Prognosis Model operates on a sample-based set of indi- 
vidual trees. Because the Budworm Model estimates defoliation on tree 
classes, a component of the damage model is required to disaggregate class- 
average defoliation estimates onto the individual trees that make up the 
classes. 



MODEL OVERVIEW 

The Prognosis Model is a distance -in dependent individual tree model, 
A user enters a sample-based list of tree measurements, data describing 
the site conditions, and program control commands. The Prognosis Model 
reads the initial data, computes some calibration statistics, and prepares 
to enter into a growth prediction loop. This loop contains the calculations 
for predicting the change for trees in the inventory list, for a given growth 
cycle. The cycle length is usually 10 years but can be changed by the user, 
Within this loop, or cycle, the following steps are followed: First, any tree 



harvesting that is desired by the user is simulated. Second, the normal 
growth rates on trees are predicted for the period of the cycle, usually 10 
years. Third, the normal mortality is estimated. At this point, normal es- 
timates of growth and mortality are known by the model. If a budworm 
outbreak is being simulated, the Budworm Dynamics Model is called. 

The calling sequence for the Budworm Model includes three major steps. 
First, the foliage biomass on individual trees is estimated. The program 
classifies these data into three tree-size classes for each host species. 
These tree classes are called Budworm Model trees to distinguish them 
from Prognosis Model sample trees. The foliage is further classified into 
three crown thirds: top, middle, and bottom. The three tree-size classes 
and three crown thirds per tree define nine foliage cells (table 2). The foli- 
age in each cell is divided into four foliage age classes (also known as foli- 
age cohorts): new, 1 year old, 2 years old, and remaining. 

The second step includes the Budworm Model's loops. The Budworm 
Model uses 1-year-long time steps; because there are 10 years in a typical 
Prognosis Model growth cycle, there are usually ten, 1-year-long Budworm 
Model cycles within a Prognosis Model cycle. During each year, the model 
predicts feeding done by budworms. The feeding reduces the amount of 
new foliage and, sometimes, older foliage as well. The foliage dynamics 
model ages the new foliage into 1-year-old foliage, and the 2-year-old foli- 
age is added to the remaining foliage. A foliage death model predicts how 
much of the remaining foliage dies each year. A foliage birth model pre- 
dicts, how much new foliage is produced the following year. The process of 
foliage birth, feeding by budworm, aging, and natural death is simulated 
for each year of the 10 Budworm Model cycles. 

The third step of the Budworm Model is to translate the budworm -caused 
damage done to the Budworm Model trees into estimates of reduced growth, 
increased mortality, and topkilling of the Prognosis Model sample trees. 
At this point, the Budworm Model logic ends and the Prognosis Model 
picks up where it left off. Other than the fact that tree growth and mortal- 
ity rates have been changed, and some trees have been topkilled, the Prog- 
nosis Model is unaware that the Budworm Model was involved. The fourth 
step of the Prognosis Model is to add increments of growth; diameter 
growth is added to diameter and height growth to height. Changes in 
crown ratio are also computed. The Prognosis Model stores the number 
of trees per acre represented by each sample tree; mortality is subtracted 
from this value. The fifth step is to summarize stand values and write re- 
ports. The process is repeated for each cycle in a Prognosis Model run. 

Sheehan and others (1989) provide another summary of the linked 
submodels and provide figures illustrating the Budworm Modeling System 
organization. The description in this paper is a little different than the de- 
scription found in Sheehan and others (1989) because this description em- 
phasizes the role of the foliage and damage models. Crookston and others 
(1990, p. 3) provide an organizational chart of the Budworm Modeling Sys- 
tem components that provides more data on the spatial and temporal scope 
of model components and their relationships to each other. 

In the following presentation, several variable names are introduced. 
Table 3 provides a summary of the variables, subscripts, and notation used 
in this paper. 



Table 3 Variables, subscripts, and notation used in this paper 



Variable 

notation 



or 



Description 



Summation over all n elements 

Foliage age class subscript, 1 = new, 2 = 1 -year-old, 3 2-year-old, and 
4= remaining 

The amount of biomass that will be aged from a younger age class to an older one 

Adjusted potential biomass 

Biomass after feeding 

Biomass before feeding 

Prognosis Model sample tree biomass computed from Moeur (1981) 

Crown foliage cell subscript deffned in table 2 

Crown length in feet 

Cumulative (5-year) defoliation 

Crown ratio 

Diameter at breast height in inches 

Periodic (10 years) growth In squared diameter (inches) 

Defoliation proportion; the exact computation depends on the use of this variable 

The function Illustrated in figure 1 

A function that predicts proportion defoliation by host, crown cell, and foliage age 
class; this function represents the Budworm Dynamics Model 

The function illustrated in figure 2 

Host subscript 

Height in feet 

Height growth In feet 

Prognosis Model sample tree subscript 

Tree-size class subscript, 1 = small, 2 medium, and 3 large 

An indexing function that returns tree size 1 for crown cells 1 , 2, or 3; tree size 
2 for crown cells 4, 5, or 6; and tree size 3 for crown cells 7, 8, and 9 

Time subscript, where 1 = initial (inventory) year within a Prognosis Model cycle. 
t usually runs from 1 to 10 

The number of years the model has run thus far. For example, if the model is In 
the fourth year of a run, T= 4 and if the model is in the tenth year T 10 

Stand density in trees per acre 

Potential biomass by host, crown cell, and age 

Proportion of expected basal area increment by host and tree size 

Proportion of expected height growth by host and tree size 

Proportion retained biomass by host, crown cell, and age 

Proportion of stem kiEled 

The number of trees per acre represented by Prognosis Model sample tree n 

Probability of topkill 

Probability of tree mortality 

The number of Budworm Model trees per unit area 



BIO n 

c 

CL n 

CMD n 

CR n 

D n 

DDS n 

DEF n 

f1 c 

f?t 

13 

h 

H n 

HG n 

n 

s 

S(c) 

t 



TPA 

PBhcat 



PRBhcat 

PRK n 

PRQB n 

PTK n 

PMRT n 



Total Foliage 
Biomass per Tree 



Dividing Foliage 
Into Crown Thirds 



The amount of total biomass on each Prognosis Model sample tree is esti- 
mated using models developed by Moeur (1981). She published several 
models that predict total foliage biomass for undamaged trees using vari- 
ables carried in the Prognosis Model. One of the models was chosen to use 
here: 

WhenD n (d.b.h. inches) < 3.5 

BlO'n = exp (B0 n + 1.468 In D n + 0.3088 In DDS n - 1.077 In H n + 

0.6908 In CL n - 0.1421 In TPA + 0.3992 In RD n ) (1) 

otherwise 

BIO' n = exp (Bl n + B2 n In CL n - 0.12975 In TPA + 0.4035 In H n ) - 

1.13178 (2) 

where 

-B/0',i = biomass is pounds per tree for tree n 

B0 m Bl n , and B2 n ~ species- specific constants for tree n (see table 4) 

D n = diameter at breast height (d.b.h.) in inches of tree n 

DDS n = periodic (10-year) increment in the square of d.b.h. for tree n 

(proportional to basal area increment) 
H n ss height in feet of tree n 
CL n = crown length in feet of tree n 
TPA =- stand density in trees per acre 
RD n ~ d.b.h. of tree n, divided by quadratic mean stand d.b.h. 

Then pounds are converted to grams; BIO n ~ 453.6 B/0' n . 

Available Data Jenson (1976) studied the foliage distribution within 
one codominant Douglas-fir tree in the Cascade Mountains of Washington. 
Table 5 presents the foliage weight for that tree by whorl and foliage age 
class. 

The whorls of Jensen's tree were arbitrarily divided into three groups 
each group corresponding to a crown third. For small trees, whorls 1 to 3 



Table 4 Species-dependent coefficients for the 
total foliage biomass equations (from 
tables 9 and 13 of Moeur 1981) 



Intercepts 


met 

D n < 3.5 
B2 


Species 1 


D n > 3.5 
BO 


D n < 3.5 
B1 


WP 


2.66607 


-1.94951 


1 .22023 


we 


3.05935 


-2.24876 


1.37600 


WL 


1.75654 


-4.73762 


1.98479 


WH 


2.65457 


-4.17456 


2.00749 


AF 


3.06017 


-1.60998 


1.32649 


LP 2 


2.622505 


-3.13488 


1 .62368 


GF 


3.11508 


-2,43200 


1 .60270 


ES 


3.30085 


-2.93508 


1.96125 


DF 


2,70587 


-2.93508 


1,25837 


PP 


2.45249 


-2.74410 


1.58171 



1 See table 1 for species list. 

Z BO for LP is a coefficient fit for whitebark pine (Pinus 

slblcaulis Engelm.). 



Table 5 Foliage biomass by whorl and four age classes for ono codominant Douglas-fir 

tree from the Washington Cascade Mountains (reorganized* data from table 9 

ofJenson 1976) 

Foliage age class 



Whorl 


New 


1-yr- 


old 


2-yr-old 


Remaining 


Total 




Pet 


2 9 


Pct 


9 


Pet 


9 


Pet 


9 


9 


1 


100 


1 




















1 


2 


85 


11 


15 


2 














13 


3 


67 


28 


28 


12 


5 


2 








42 


4 


54 


24 


38 


19 


8 


4 








49 


5 


41 


36 


34 


30 


18 


16 


7 


6 


89 


6 


39 


43 


32 


35 


19 


21 


10 


11 


109 


7 


34 


60 


32 


56 


23 


40 


11 


19 


176 


8 


31 


56 


28 


50 


23 


41 


18 


32 


181 


11 


119 


49 


22 


57 


19 


49 


40 


104 


260 


14 


110 


29 


21 


61 


18 


52 


51 


198 


290 


17 


177 


23 


16 


52 


15 


49 


62 


201 


324 


20 


200 











1 


1 


99 


21 


22 



1 Some of the percentages reported by Jenson are slightly altered so that each row sums to 100 

percent. 
2 Grams dry weight. 



represent the top third, whorls 4 to 6 represent the middle third, and 
whorls 7, 8, and 11 represent the bottom third. This represents an 11-year-old 
tree, which is a reasonable assumption for an average tree between and 
7 m tall. For medium and large trees, whorls are assigned to thirds as fol- 
lows: top third is whorls 1 to 5, middle third is whorls 6 to 11, and bottom 
third is whorls 14 to 20, Table 6 illustrates the percentages of the total 
foliage in each crown third of Jenson's tree using these arbitrary whorl 
classifications. 

Schmid and Morton (1981) reported mean percentage foliage area and 
biomass in various crown levels of Douglas-fir and white fir. They sampled 
eight open-grown Douglas-fir and 10 white fir from Capulin Mesa near Los 
Alamos, NM. Garolin and Coulter (1972) list percentages of foliage area for 
crown thirds of Douglas-fir trees sampled in Oregon, All these data, in- 
cluding the data from table 6, are displayed in table 7. 

Using a Volume Approach Another method for computing the per- 
centages of total foliage by thirds is to classify a tree into a crown shape, 
compute the volume of each third, and express the volume as a percentage 
of the total. It is useful to compare the percentages presented in table 7 to 
those derived analytically for a given crown shape. A convenient crown 
shape for this comparison comprises of a cylinder of radius r and height 2h 
for the bottom two-thirds, and a cone of radius r and height h for the top 
third. The crown length of our hypothetical tree is 3/1. Using this method, 
the crown volume in the top is approximately 14 percent, for the middle 
and bottom thirds it is 43 percent. 

Technique Used There is a notably strong similarity between all of 
the estimates for medium to large Douglas-fir trees. The similarity be- 
tween percentage of foliage area and percentage foliage weight is strong 
in Schmid and Morton's (1981) data. This lends support to the idea that 
measurements of percentage foliage area by Carolin and Coulter (1972) can 
be considered as additional data. The percentages derived from Jenson's 



Table 6 Biomass by crown third for Jenson's (1976) tree using arbitrary whorl 

classification to define crown thirds 



Tree 

class 


Crown third 


Top 


Middle 


Bottom 


Total 




Pet 


a 


Pet 


9 


Pet 


9 


Pet 


9 


Small 


6 


56 


27 


245 


67 


617 


100 


918 


Medium and 


















large 


12 


191 


47 


723 


41 


637 


100 


1,151 



Table 7 Summary of data on percentage of foliage by crown thirds 



Crown class 


Description and source 


Top 


Middle 


Bottom 


Small DF trees, percentage of biomass 








(table 6) 


6 


27 


67 


Medium and large DF trees, percentage 








of biomass (table 6) 


12 


47 


41 


Medium DF trees, percentage of biomass 








(Schmldand Morton 1981) 


16 


55 


29 


Medium DF trees, percentage of foliage 








area (Schmid and Morton 1981) 


15 


52 


33 


Medium WF, percentage of biomass 








(Schmid and Morton 1981) 


24 


49 


27 


Medium WF, percentage of foliage area 








(Schmid and Morton 1981) 


24 


44 


32 


Medium DF, percentage of foliage area 








(Carolin and Coulter 1972) 


12 


36 


52 


Volume approach described in 








this paper 


14 


43 


43 


Values used In the model 








Small trees 


5 


30 


65 


Medium and large trees 


15 


45 


40 



Dividing Foliage 
Into Age Classes 



(1976) data, the volume approach, and the direct measurements are similar. 
The choice of arbitrary divisions in the whorl numbers of Jenson's tree yielded 
percentages that are not grossly different from the other measurements. 

I believe that errors in estimating proportions of foliage in crown thirds 
are relatively unimportant to meeting the objectives of this foliage model. 
Therefore, I computed a simple average of the percentage data and ad- 
justed the numbers to even 5 percentage points such that their sum equals 
100 percent (see the bottom of table 7). 

Within each crown third, the percentages of the foliage that are new, 
1-year-old, 2 years old, and the remaining need to be computed. 

Available Data Jenson's (1976) DF (see table 1 for species abbrevia- 
tions) tree data (table 5) can be organized into the whorl groups used to di- 
vide the foliage into crown thirds. The relative amount of foliage by foliage 
age class can then be computed for the crown thirds. These data are dis- 
played in table 8. 



Table 8 Biornass by foliage age class for each crown third derived from Jensen's (1 976) data, 
table 5 



Tree 
size 


Whorl 


Crown 


Foliage age class 


class 


numbers 


third 


New 


1-yr-old 


2-yr-old 


Remaining 








Per 


9 


Pet 


9 


Pet 


9 


Pet 


9 


Small 


1-3 


Upper 


71 


40 


25 


14 


4 


2 










4-6 


Middle 


42 


103 


34 


84 


17 


41 


7 


1/ 




7-11 


Bottom 


27 


165 


26 


163 


21 


130 


25 


155 


Medium, 


1-5 


Upper 


52 


100 


33 


63 


12 


22 


3 


6 


Large 


6-11 


Middle 


29 


208 


27 


198 


21 


151 


23 


166 




14-20 


Bottom 


8 


52 


18 


113 


16 


102 


58 


370 



Computing Initial 
Biomass 



Hatch and Mika (1978) reported the percentage of new foliage on 18-inch- 
long midcrown sample branches from DF and GF. For DF, the mean was 
27 percent and the standard deviation was 13 percent. For GF, the mean 
was 35 percent and standard deviation was 7 percent. These data do not 
include an estimate of the percentage of new foliage on the crown third. 
The 18-inch-long samples may or may not include all needle cohorts in the 
midcrowns. However, the mean Douglas-fir estimate from Hatch and Mika 
is 27 percent, and the percentage for Jensen's tree is 29 percent. This fact 
supports the data derived from Jenson's tree. 

No other data on foliage profiles were found for the host species. 

Technique Used Table 9 displays the percentages used in the foliage 
model to partition foliage in each crown third. These values are not identi- 
cal to those in table 8; they are slightly different due to arbitrary rounding 
to ensure the two totals sum to 100 percent and, in the case of upper 
crowns, a small amount of remaining foliage is present. The reason the re- 
maining foliage was added to the upper crowns of small trees is related to 
computational convenience within the dynamic part of the foliage model 
(described below). 

Potential Biomass Recall that Moeur's (1981) models (equations 
1 and 2) estimate total foliage biomass on an undamaged tree, and that 
those estimates are converted to grams in variable BIO n > LetPJ5ft ca be the 
total biomass per tree for host h, foliage cell c (see the nine foliage cell 
codes defined in table 2), and foliage age class a (1 = new, 2 = 1-year-old, 
3 = 2-year-old, and 4 = remaining foliage). Then; 



PC C PROB n )/Z n PROB n 



(3) 



where 



= denotes a summation over sample trees of host h 

~ proportion of foliage of age class a within foliage cell c (see 

table 9) 

= proportion of the tree biomass in crown class c (see table 7) 
= trees per acre represented by tree n. 

Adjusted Potential Biomass In a budworm-infested forest, the actual 
biomass differs from potential biomass because budworm feeding alters the 
tree's foliage profile. The amount of old foliage can be greater than, normal 
due to trees' tendency to retain old needles longer than normal during 



"Aca 

PC C 

PROBn 



Table 9 Percentage of new, 1 -year-old, 2-year-old, and remaining 
foliage by crown third for undamaged trees 



Tree 

size 
class 


Crown 
third 




Foliage 


age class 




New 


1-yr-oid 


2-yr-old 


Remaining 


Small 


Upper 
Middle 


70 
45 


24 

30 


4 
20 


2 
5 




Bottom 


35 


30 


25 


10 


Medium 


Upper 

Middle 


50 

30 


35 
25 


12 
20 


3 
25 




Bottom 


10 


20 


20 


50 


Large 


Upper 
Middle 


50 
30 


30 
25 


15 
20 


5 
25 




Bottom 


10 


20 


20 


50 



periods of defoliation. A. Van Sickle (in Thomson 1979, p. 12) reported this 
phenomenon. 

The amount of new foliage produced by a tree is partly a function of the 
defoliation during the preceding year. Defoliation- caused reductions in 
food reserves reduce the number of needles in buds and the weight of new 
foliage in balsam fir (Mayfield 1984). Thomson (1979) hypothesized that 
the same phenomenon occurs in Douglas-fir and suggests that new foliage 
depends partly on defoliation of new foliage in the previous year. Thomson's 
suggested function was parameterized for the tree classes used in the Bud- 
worm Model (fig. 1). The parameters were chosen in 1982 to match a hy- 
pothesis proposed by some of the scientests listed in the acknowledgments 
of this paper. 

The adjusted potential biomass before feeding (APBh cat } for time t - 1 
(the inventory year) is computed from PB hca , the functions illustrated in 
figure 1, and inventory- time estimates of previous defoliation; 



APB hcat = PB hca fl c (PRB hc(a+1)i ) 
where 



(4) 



t = proportion retained biomass for host h, crown c, foliage 

age class a, at time t 
h, - host species 

c = foliage cell code (table 2) 

a = foliage age classes 1, 2, and 3 

t =1 (inventory time) 

fl c ~ function illustrated in figure 1, indexed by foliage cell. 

The inventory-time estimates of the proportion of retained biomass are 
entered by the model user (Crookston and others 1990, p. 7). For example, 
a value of 0.80 for a given age class in a given crown of a given host is an 
indication that 20 percent of the biomass on that branch was missing (pre- 
sumably eaten or destroyed by budworms). Values of PRB for the remain- 
ing foliage class may be greater than 1.0 indicating that more than normal 
levels of remaining biomass were present at inventory time. 

The adjusted potential biomass corresponds to the amount of foliage an 
observer would measure on a tree that may have been defoliated in the re- 
cent past but is not currently defoliated. Within the Budworm Model, APS 
is the denominator used to compute percentage defoliation. 
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Proportion of potential 

new foliage produced 

1.0 



0.8- 



0.6- 




Tree size and 
crown level 



Small 





Medium 




Large 



0.2 0.4 0.6 0.8 

Proportion of retained biomass 



Figure 1 Relationship between proportion of retained biomass and 
the proportion of the potential new foliage a tree can produce. Line A 
applies to all crown levels of small trees, line B applies to the tops of 
medium and large trees, tine C applies to the middle of medium and 
large trees, and line D applies to the bottom of medium and large trees. 



1.0 



Biomass Before 
and After Feeding 



Annual Foliage 
Dynamics 



Biomass before feeding (BBF) is the amount of foliage available to bud- 
worms. The initial computation of BBF is 

BBF hcal = APB hcat PRB hcat (5) 

for all four age classes; the other variables are as defined above. The Bud- 
worm Model contains component models that convert the new foliage class 
into bud numbers and simulates budworms eating buds developing new fo- 
liage. Spring foliage development rates are predicted using models devel- 
oped by Beckwith and Kemp (1984). Biomass after feeding (BAF) is repre- 
sented by 

BAF hcat = BBF hcat f2 hcat (6) 

, where f2 represents proportion defoliated by the Budworm Model on host A, 
crown c, age a, and in year t. There is an option in the Budworm Model 
(Crookston and others 1990, p. 7) that allows users to circumvent the Bud- 
worm Dynamics Model and directly enter values for f2. This option was 
used in preparing examples of behavior for this paper because it allows 
users to isolate the behavior of the foliage dynamics and damage compo- 
nents from the behavior of the Budworm Dynamics Model. 

Foliage estimates for time t - 2, 3,. ..,10 are computed as described below, 

Recompute Proportion of Retained FoliageAfter a year of feeding, 
a new estimate ofPRB can be computed: 

PRB hcat = BAF hcat /APB hcat (7) 
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Compute Remaining Foliage The amount of remaining foliage de- 
pends on the amount being added from the 2-year-old foliage class as this 
foliage ages, and the amount that dies. The idea of retaining old biomass 
described earlier is incorporated in this model. 

Part I: The foliage model assumes that whatever the current status of 
remaining foliage, the tree tends toward normal levels when defoliation 
stops. When PRBfait is greater than or equal to 0.80 (new foliage defolia- 
tion is less than 20 percent), the model assumes that defoliation of new fo- 
liage during year t has essentially stopped. The amount of remaining foli- 
age present when defoliation stops can be greater, equal, or less than the 
amount found on undamaged trees. Excess old foliage (40 percent per 
year) dies quickly; below-normal amounts are replaced by aging 2-year-old 
foliage. 

The aging of 2-year-old foliage into the remaining foliage age class is 
complicated by the differences in the potential biomass in the 2-year-old 
and remaining class. The proportionate difference between these classes is 
used to scale the amount of biomass that moves from the 2-year-old to re- 
maining. LeiAGEkc(a+i)t be the amount of foliage aging from age class 3 to 
class 4, in host h, crown c, at time t: 

AGE hc4t = BAF hc3l (PB hc4 /PB hc3 ) (8) 

The rest of Part I of the remaining biomass calculations are as follows : 
IfBAF hc4t >PB hc4 ihen 

BBF hc4(t+1) = PBhc4 when (BAF h c4t - PB hc4 }/PB kc4 < 0,10, otherwise 

BAF hc4t - (0,4 (BAF hc4i -PB hc4 )) (9) 

If BAF hc4t < PB hc4 , let X BAF hc4t + AGE hc4t then 
BBFh c4(t+ 2) = X when X < PBhct and PBh C 4 otherwise (10) 

Part II; When PRB hclt is less than or equal to 0,80 (defoliation is over 
20 percent), the amount of remaining foliage is driven by the tree's ten- 
dency to reduce foliage mortality, 

Three mortality effects are hypothesized. First, if the remaining biomass 
is greater than or equal to the potential biomass summed over all foliage 
age classes, as much biomass dies each year as is added from the 2-year- 
old class. This means that a tree cannot retain more biomass in the re- 
maining category that it could potentially hold over all foliage age classes. 
This becomes an upper bound on the remaining foliage class. Any foliage 
accumulating beyond this point is killed by the model. 

Second, when a tree is not being defoliated, the amount of old foliage dy- 
ing approximately equals the amount being added to the old (remaining) 
foliage class. Thus, the intrinsic mortality of the remaining foliage would 
be equal to the amount of 2-year-old biomass that is aging into the remain- 
ing class. In this model, the intrinsic mortality is reduced as new foliage 
defoliation increases (fig. 2). The rate is further adjusted to account for the 
fact that the amount of old foliage dying is much smaller than the amount 
being added to the remaining foliage class in small trees. The rates used 
are specific to the host species. As modeled, Douglas-fir and spruce have a 
stronger ability to retain old foliage when new foliage is lost than grand fir, 
white fir, and subalpine fir. 

The third mortality effect is directly related to total tree defoliation. It 
accounts for the idea that a tree's ability to retain old needles is hampered 
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0.0 



0.2 0.5 0.0 

Defoliation of new foliage 



Figure 2 Relationship between new foliage 
defoliation and the proportion of the foliage 
entering the old foliage class that is retained. 
Old needle retention is stronger in Douglas-fir 
and spruce than In the firs. 



due to branch and needle drying and sun scald coincident with opening up 
the canopy to additional light. Evidence that this might happen is limited 
to the fact that when the model was constructed without this effect, surviv- 
ing heavily defoliated trees retained old needles indefinitely unless the old 
needles were actually eaten by budworms. 

This last effect is modeled by monitoring total tree defoliation. When to- 
tal defoliation exceeds 40 percent, 15 percent of the remaining foliage dies. 

A formal statement of the model follows; 

LetXA/t = 1.0 if total tree defoliation in year t is less than 40 percent and 
0,85 otherwise. Let 

XAB hc4 t - (BAF kc4l +AGE hc4 tf3(PRB hcl ))XM t (11) 

where f3 is the function displayed in figure 2.. 
Finally, 

BBF hc4 (M) = XAB hc 4t when XAB hc4t < la PB hca , and la PB hca (12) 

otherwise. 

Compute 1- and 2-Year-Old Foliage The 1- and 2-year-old foliage is 
computed by aging the younger foliage. The process moves the 1-year-old 
to the 2-year-old then the new to the 1-year-old. For a = 2 and 3 



BAF hc(a ~ 1)t (PBhJPB^.!)} (13) 

Compute New Foliage The amount of new foliage each year is mod- 
eled as a function of the potential new biomass an undamaged tree can 
grow and the defoliation the tree has undergone, This is the same model 
as that used for initial conditions except that the measure of defoliation ia 
a weighted average of two measurements. The first measurement is the re- 
tained biomass after feeding divided by biomass before feeding. The second 
measure of defoliation is the biomass after feeding divided by the potential 
biomass expected on a healthy tree. 
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Neither measurement has been directly related to new foliage develop- 
ment in field studies. However, both have some theoretical basis. The first 
measure is given a weight of 0.40 and the second measure a weight of 0.60. 
The second weight is greater than the first, implying that second measure 
more correctly measures the ability of the tree to regenerate foliage. The 
following function is used to compute new foliage in the next period: 

BBFkcXt+v = PBhd /2(0.4 PRB hc2 + 0.6 (BAF hc2t /PB hc2 )) (14) 

Periodic Foliage Dynamics The annual foliage dynamics described 
above are used to model year-to-year dynamics for approximately 4- to 10- 
year-long periods. At the end of the period, the model computes a final es- 
timate of proportion retained biomass and stores it for use at the beginning 
of the next period to recompute initial conditions. 

The Prognosis Model estimates tree growth and mortality at every pe- 
riod. Therefore, at the beginning of the next period, some of the trees that 
were in the small- or medium- tree -size class will grow into the next class. 
For example, let's say your stand has a few small trees in it. The medium 
class is empty. At a period boundary, the trees might "move" to the next 
larger class. The foliage model has a provision in it that traces the move- 
ment of defoliated trees from one class to another and properly adjusts the 
estimates of proportion retained biomass to reflect the movement. 

Foliage Model The behavior of this model is examined in two respects. The first con- 

Behavior cerns the changes in several state variables under different defoliation pro- 

files, The period will be held at 10 years, and the annual behavior will be 
exercised, illustrated, and discussed. The second concerns the periodic be- 
havior. In this model, the period length is user-input, and the structure of 
the model differs greatly on period boundaries versus annual time steps. It 
is important, therefore, to understand the impact changing the period 
length has on model behavior. The same stand and defoliation profiles are 
used to examine both the annual and periodic behavior. 

Description of Example Site-The stand data used in the examples of 
behavior were collected from Bear Gulch, about 20 km south of John Day 

TT Q ^ 6 u Slte 1S a bng and Steep north - facm g slope visible to travelers of ' 
U.fa. Highway 395. The south-facing slope to the north of the study site 
contains little timber; the same is true for the south-facing slope located 
south of the study site. Therefore, the site is bounded on the north and 
south by a significant area of essentially nontimber lands. The area is a 
Douglas-fir site that contained many large ponderosa pine trees until they 
were logged in the mid-1980's. The remaining stand of mostly Douglas-fir 
was heavily defoliated by western spruce budworm soon after the site was 
logged Table 10 lasts several characteristics of the example site (this site 
is used in Crookston [1991] as well). 

Three defoliation profiles are applied to the stand (table 11). For case 1 
95 percent defoliation of the new foliage and 40 percent of the 1-year-old fo- 
liage is applied for each year of a 10-year simulation. Case 2 is the same 
except that defoliation is stopped after 5 years. Case 3 uses variable defo- 
liation which more nearly typifies actual experience. All tree-size classes 
and foliage cells are treated identically. 



t re plotted over time to us 

the behavior of the foliage dynamics models. All of the figures use the 

same format. The variable labels (figs. 3-6) and definitions are as follows* 
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Table 10 Mensurational characteristics of 
the Bear Gulch study site 

Attribute 1 Value 

Number of trees 

per hectare 1,003 

per acre 406 

Basal area 

m 2 /ha 25.0 

ft 2 /acre 109 

Top height 

meters 23.2 

feet 76.3 

Merchantable volume 

m 3 /ha 157.9 

ft 3 /acre 2,256 

Percent Douglas-fir 89 

Percent ponderosa pine 1 1 

1 Conversion factors: 1 ft = 0.3048 m, 
1 inch a 2,54 cm, 1 acre = 0.4047 ha, 
1 ft a /acre basal area = 0.2296 m 2 /ha, 
1 ft 3 /acre merchantable volume = 0.0700 m 3 /ha. 



Table 11 Percentage defoliation values used in the three test cases 
Case 1 Case 2 Case 3 



Year 


New 


1-yr-o|d 


New 


1-yr-old 


New 


1-yr-old 


1 


95 


40 


95 


40 








2 


95 


40 


95 


40 


100 


10 


3 


95 


40 


95 


40 


20 





4 


95 


40 


95 


40 


90 


20 


5 


95 


40 


95 


40 


20 





6 


95 


40 








95 


10 


7 


95 


40 








50 





8 


95 


40 














9 


95 


40 














10 


95 


40 















Adj-Pot Adjusted potential biomass in grams per tree on an average 
tree prior to defoliation each year. Let WT hs be the number 
of trees per unit of area of host h and tree size s and S(c) be a 
function that provides a tree size s given a crown index c. 
Then; 



Adj-Pott ^dhU^APBh^WThs^/Ihlc WT hS(c ) (15) 

Tot-Post Adj-Pot after feeding. 

Rem-Post Remaining biomass in grams per tree after feeding on an 
average tree. 

Rem~Post t = (S*Sc BAF hc4t WnsiyVMe WT hS ( c ) (16) 

New-Pre New biomass in grams per tree before feeding on an average 
tree. 

New-Pre t (S ft lc BBF hclt WT hs(o ))/IhIc W?hS(c) (17) 

New-Post New-Pre after feeding, 
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Year 

Figure 3 Foliage dynamics for case 1 (table 1 1 ) 
indicated by AdJ-Pot, the adjusted potential 
blomass on an average tree; Tot-Post, which Is 
Adj-Pot after defoliation; Rem-Post, the blomass 
of remaining foliage after defoliation; New-Pre, the 
new biomass before defoliation; and New-Post, 
the new biomass after defoliation (refer to the text 
for complete definitions). 




New-Post 



Figure A Foliage dynamics for case 2 
(table 11) illustrated by plotting five 
indicator variables over time (see fig. 3 
and text). 
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9 10 

Year 

Figure 5 Foliage dynamics for case 3 (table 1 1 ) 
3d by plotting five indicator variables over 

io fin Q and imrtV 



./ ./Tot-Post 
. ././RenvPost 
/ ./Now-Pre 

^-/New-Post 
m 



Figure 5 Foliage dynam 
illustrated by plotting five I 
time (see fig. 3 and text). 




Figure 6 Foliage dynamics for case 2 run in one 
10-year-long Prognosis Model cycle (APot-10, 
TPost-1 0, and NPre-1 0) and in two 5-year-long 
Prognosis Model cycles (APot-5, TPost-5, and 
NPre-5). The difference between pairs of lines is 
due to changes caused by the Prognosis Model 
cycle boundary. Label associations with those 
used in figures 3-5 are: APot = Adj-Pot, TPost = 
Tot-Post, NPre New-Pre. 



17 



Annual Behavior Annual behavior is examined by running the model 
for 10 years without exercising the periodic component of the model. Figure 3 
corresponds to case 1 (table 11), figure 2 to case 2, and figure 4 to case 3 
The initial conditions (year 1, figure 3) indicate that there are 5,200 g per 
tree of adjusted potential biomass, and because there had been no previous 
defoliation, the adjusted potential (APB) and the potential (PB) are the 
same. There are just over 1,000 g per tree of new foliage available for 
the first generation of budworms to feed on. 

( In case 1 (fig. 3) the adjusted potential and remaining foliage age class 
rises for the first year following heavy defoliation. This rise is due to the 
model component that increases old needle retention at the onset of defo- 
liation. As high defoliation continues, the needles from the remaining class 
cue, resulting m a steady decline in foliar biomass in the remaining class as 
well as on the entire tree. New foliage prior to defoliation is over 1,000 g in 
the first year then drops to below 300 g in subsequent years. This is due to 
the model components that represent the effect defoliation has on new foli- 
age production. 

In case 2 (fig. 4), the same trends exhibited in figure 3 are evident for the 
farstS years. After the fifth year, the foliage recovery is displayed. New fo- 
liage biomass returns to preoutbreak levels in about 3 years. Foliage in 
the remaining class returns to preoutbreak levels quickly because the lev- 
els were not far from normal when the defoliation stopped. All of the indi- 
cators return to normal before the end of the 10-year cycle 

Case 3 (fig 5} illustrates the interaction of the foliage model components. 
The pulses of defoliation caused the foliage model to respond by reducing 
the amount of new foliage in the next year; remaining needles were re- 
tained and the total in this foliage age class increased. All of the effects 
of defoliation were followed by the effects of recovery when the defoliation 
levels were low. 

Periodic Behavior-All three cases were run using one 10-year-long 
Prognosis Model cycle and then using two 5-year-long cycles. Figure 6 dis- 
plays the results of running case 2. Three of the five indicator variables 
are illustrated as pairs of lines in figure 6. All three are closely related 
suggesting that the foliage dynamics model behaves similarly when there 
is a cycle boundary in the middle of an "outbreak" as compared to running 
one 10-year-long cycle. The other cases are not illustrated because the 
cycle boundary effect was even smaller in those cases than in case 2 
Ihe effect of tree growth between cycles is minimal during an infestation. 
I his is partly due to the fact that defoliated trees do not grow well In a 
model run of lightly defoliated trees that are fast growing, a discontinuity 
in the indicator variables caused by tree growth will be more evident Fig- 
ure 6 shows a slight difference in the lines labeled TPost-5 and TPost-10 
starting at year 5. This difference is due to tree growth. 

The effects of a cycle boundary in the middle of a Budworm Model run 
are not strong as far as the foliage dynamics model is concerned There 
may be other subcomponents of the Budworm Modeling System that are 
less well behaved in this respect. 

THE TREE DAMAGE MODEL 

Defoliation of trees reduces diameter and height growth, increases mor- 
tality rates, causes topkilling, and influences the success of regeneration 
establishment. More is known about the effects of defoliation on tree 
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Using the Beta 
Distribution 



growth than is known about the foliage dynamics of defoliated trees. How- 
ever, each study that relates the effects of defoliation on tree growth and 
development uses a different measurement of defoliation. Therefore, part 
of the damage model is devoted to computing defoliation metrics. 

The Budworm Model represents trees as classes: small, medium, and 
large; and the Prognosis Model represents trees as individual sample trees. 
The defoliation predicted on a given Budworm Model tree class ia consid- 
ered an average for all of the Prognosis Model trees in the class. In natural 
stand conditions, all trees within a tree-size class do not have identical de- 
foliation patterns (McDonald 1981). To represent this variation, the simu- 
lated average defoliation for each tree class is distributed to the Prognosis 
Model sample trees that make up the class following the Beta distribution. 
The reasons for using the Beta distribution and the way the distribution 
parameters are computed are described next. 

Data To develop growth equations for budworm-defoHated small trees, 
Ferguson (1988) measured the defoliation in the tops, middles, and bottoms 
of individual trees on each sample study site. Ferguson's data, augmented 
with a small portion of data taken in Oregon, were used to parameterize 
the Beta distribution used in the damage model. The trees on each study 
site were classified into the tree-size classes used by the Budworm Model 
(table 2). Figure 7 plots the variance of the individual tree defoliation 
measures over the mean for each foliage cell and study site. Cases with 
less than four observations were deleted. 



0.00- 



0.0 



T M 




0.2 0.4 0.6 O.f 

Class average defoliation 



1.0 



Figure 7 The class variance of individual tree defoliation 
measures plotted over the class average. Symbols are T top 
crown third, M = middle third, B = bottom third, and A average 
of T, M, and B. The solid line represents the function used by 
the model to compute a class variance given the class average. 



19 



w = 1 - 



As the mean defoliation increases toward 50 percent, the variance also 
increases, and when the mean exceeds 50 percent, the variance diminishes 
toward zero. Theoretically, the variance of proportion data must be zero 
when the mean is zero or 100 percent. The peak variance could be shifted 
to the right or left of 50 percent (the underlying distribution need not be bi- 
nomial). These data indicate the peak is near 50 percent and that the dis- 
tribution is not asymmetrical. (The line in fig. 7 is discussed below.) 

Justifications The justifications for using the Beta are as follows. 
First, 100 percent of the probability density function (p.d.f.) is contained on 
the interval (0-1). Second, the shapes of the p.d.f. are appropriate for these 
data. The shape parameters (v and w) for the p.d.f. are estimated from the 
sample mean, x and variance, s 2 (Hastings and Peacock 1975): 

(18) 
1) (19) 

Figure 8 illustrates several shapes of the Beta where the mean is set to 
0.5 (50 percent) and the variance is varied from 0.02 to 0.12. When the 
variance is low, most of the probability density is near the mean. When 
the variance is high, the probability density is pushed toward the tails. 
For defoliation data, the Beta appropriately implies that if a class of trees 
has a mean proportion defoliated of 0.5 (50 percent) and a variance of 0. 12, 
then the number of trees defoliated between and 0.1 (10 percent) plus 
those defoliated between 0.9 and 1.0 would exceed the number defoliated 
between 0.4 and 0.6. If the variance is 0.02, the number defoliated be- 
tween 0.4 and 0.6 would exceed the number found in the tails of the 
distribution. 

Calculating the Parameters The data in figure 7 clearly indicate 
that the class variance is much lower when the mean is low or high than 
it is when the mean is 0.5. The line in figure 7 illustrates a reasonable, al- 
beit very conservative, represention of the relationship between class vari- 
ance and class mean. This function is taken from the Weibull p.d.f.; 




Fractlle x 

Figure 8 The Beta probability density function 
for several variances with the mean Is held 
constant at 0.5. 
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s 2 = a V(cx c ~ 2 exp (~x c )) (20) 

where a = 0.9304, c = 2.7, Vis the maximum variance, and x is the sample 
mean. This function is illustrated for V = 0.06 and 0.03 in figure 9 (the rea- 
son V = 0.03 is displayed will become clear in following sections). Although 
this function is nearly symmetrical when c = 2,7, 1-JE is used when x is 
greater than 0.5 to guarantee symmetry. 

Figure 10 illustrates shapes of the Beta for various means whereby the vari- 
ance is computed as a function of the mean using equation 20 where V= 0.06, 
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Figure 9 The class variance as a function of the class 
average for two values of maximum variance. 
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Figure 10 The Beta probability density function for several means 
whereby the variance is computed using the function displayed In 
figure 9, with maximum variance set to 0.06. For means 0.01 , 
0.05, and 0.1, only part of the defined Interval of the fractlle Is 
shown. Complementary curves exist but are not illustrated. For 
example, the curve for mean = 0.8 is the mirror image of the one 
shown for mean = 0.2. 
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In the Budworm Modeling System, a defoliation value for an individual 
tree within a given tree class is generated from the Beta distribution with 
shape parameters derived from the relationships used to plot figure 9. 
Generation of Beta deviates given real shape parameters is done using 
the generator of Cheng (1978). In the following text, the notation 

x = P<,V) (21) 

is used to denote that A; is a Beta random variable with shape parameters 
computed from the mean and maximum variance (equations 18-20). 

Other than described above, the Budworm Modeling System does not ex- 
plicitly represent variation between individual trees. The Budworm Model 
trees are averages representing a range of genetic and phenologic charac- 
teristics between trees in each model class. The approach taken in the 
damage model represents the empirical evidence that all trees within a 
class are not alike with respect to defoliation. In choosing a conservative 
value for the maximum variance, the model is programmed to represent 
only a portion of the possible variance as seen in figure 7. If more informa- 
tion regarding the between-tree variability were available, a more complete 
representation of this variance may be justified. 

Some of the damage models are functions of sums or averages of several 
years of defoliation. Ferguson's (1988) data had repeated measures, as did 
the data from Oregon. A similar scatter diagram to that displayed for fig- 
ure 7 for the average of several years' defoliation, showed that the maxi- 
mum class variance of the 4- to 5-year average is approximately half the 
variance of the individual year values. Therefore, a maximum variance of 
0.03 is used where multiyear averages are used. 

Topkill Topkill occurs when live tree height is reduced. The damage model rep- 

resents topkilling in two phases. First, the probability of topkilling is com- 
puted. This probability is compared to a uniform random number (0, 1), 
and if the probability of topkill is greater than the random number, the 
tree is topkilled. Second, the amount of topkill is computed. The models 
used in this two-step process are independently derived and applied. The 
final component of the topkill logic concerns determining if topkilling per- 
manently damages a tree. That is, the model assumes that trees with tops 
killed down to a 4-inch diameter will never produce merchantable volume 
above the point of topkill. If the diameter is less than 4 inches, the tree 
may recover a merchantable top. 

Small Trees Ferguson (1988, p. 6) published a "dieback" model for 
small trees. The probability of topkill model (fig, 11) is a function of defo- 
liation, tree height, and crown ratio. His published model includes two de- 
foliation terms with unequal constants. In the Budworm Model, these 
terms are combined to produce the following; 

PTK n = 1/(1 + exp( - (-2.75817 - 0.027635 CR n + 3.709 DEFg* + 

0.0488 #))) (22) 

where 

P7!BT n = probability of topkill for Prognosis Model sample tree n 
CR n = crown ratio (percent) of tree n 
H n - height in feet of tree n 
DEF n - defoliation of tree n. 
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Figure 11 Probability of small-tree topkill plotted over 
defoliation proportion for three tree heights (modified from 
Ferguson 1988). 



In this case, DEF n is defined as follows. A running average proportion of 
new-foliage-retained biomass in the tops of trees is computed by host and 
tree-size class during each year of the model run (AW^r, where Tis the 
number of years the model has run thus far). The running average requires 
that the Budwonn Model has run at least 3 years. The minimum value of this 
running average is stored. Its complement is the maximum average defolia- 
tion over the length of the Prognosis Model cycle. Ferguson's original model 
used average defoliation over the entire tree rather than average over just the 
tops of trees. The data used to fit the model were composed mostly of small 
trees where the defoliation is similar over all crown thirds. Because the model 
is applied to larger trees that experience differential defoliation between crown 
thirds, defoliation at the tops of trees was used in equation 22, 
Formally stated, AVT h iT is 1.0 for T = 1 and 2. For 2 < T < 6, 

AVT hlT = I t PRB hlll /min(T, 5) (23) 

Only 5 years of defoliation are counted. That is, when T > 5, the summa- 
tion Z ( starts at T- 4, rather than 1. 

At end of the Prognosis Model cycle, the min is found and then DEF n is 
computed: 



min (AVT hu> AVT hl2 ,,, t , AVT hlT ] and 



(24) 



DEF n = fK(l-Amr M ), 0.03) (25) 

where tree n is of host h and a small tree (size s - 1). DEF n is then used in 
equation 22 to compute PTK n , , 

If PTK n is greater than a uniform random number, the tree is actually top- 
killed. The amount of topkill is computed by a separate function. Ferguson 
(1988, p. 7) reported no success in building a function that relates the amount 
of topkill (he uses the term dieback) to the amount of defoliation, nor to 
any other predictor variables. Therefore, he characterized his data with 
a Weibull distribution function. His recommended model is used in the 
Budwonn Model exactly as proposed. 
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Figure 12 illustrates Ferguson's function. The way this function is ap- 
plied in the model is simple. A uniform random number between and 1 
is drawn for each sample tree. The random number is a point on the axis 
labeled proportion of trees in figure 12. The corresponding point on the 
axis labeled proportion of stem killed is found, thereby completing the com- 
putation. Formally stated, let u n be a uniform random number on (0, 1) for 
tree n; then the portion of the stem killed is 

= 0.05684 (-In u n F 036 (26) 



Medium and Large Trees Ferguson's (1988) models cover 
regeneration-sized trees, which were less than 7 m (23 ft) in height. 
For larger trees another set of models were conceived by Bruce Hostetler, 
Tommy Gregg, and the author in a 1988 meeting held at Hostetler and 
Gregg's office, Forest Pest Management, Pacific Northwest Region, Forest 
Service, Portland, OR. These models represent the professional judgment 
of Hostetler and Gregg; their equation form is a modification of the Weibull 
p.d,f. The probability of topkill is computed and compared to a uniform 
random number as done for small trees. For trees that are topkilled, the 
amount of topkill is computed using a separate model. Both models are 
functions of defoliation, computed exactly the same way as done for small 
trees (equations 23-25). Figure 13 illustrates the probability of topkill 
model: 

PTK n = Bl(l - exp (- (B2(DEF n + 1)) B3 )) (27) 

where El = 0,96, B2 = 0.65, and B3 = 14, for species GF, AF, or WF; 
Bl = 0.90, B2 = 0.60, and B3 = 11 for DF and ES. 

Figure 14 illustrates the proportion of stem killed (PRK n ) model. The 
model is exactly the same as equation 27, except for species GF, AF, and 
WF; Bl = 0,0960, B2 = 0.717, and B3 = 7.92; and for species DF and ES: 
Bl = 0.2545, B2 = 0.641, and B3 = 6.60. 
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Figure 12 Proportion of stem killed plotted over 
proportion of small trees that are topkilled (from 
Ferguson 1988). 
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Figure 13 Probability of medium- and large-tree topkill 
plotted over defoliation proportion. The line conceived 
for GF is used for AF and WF; the line conceived for DF 
is used for ES. 
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Figure 14 Proportion of medium and large tree stems killed 
plotted over defoliation proportion. The line conceived for GF 
is used for AF and WF; the line conceived for DF Is used for ES, 



Height Growth 



For trees that have not been topkilled, budworm-caused defoliation can 
result in reduced height growth. For small trees, Ferguson's (1988) height 
growth model was adapted to estimate the effects of defoliation, For me- 
dium and large trees., Nichols' (1988a, b) models are used, 

Small Trees Ferguson (1988, p. 7) published a small-tree height 
growth equation of the general form In HG = X; B/JX), where BI are regres- 
sion coefficients and^ are independent variables, The independent vari- 
ables include two measures of defoliation: ^i is for the first 2 years of the 
6-year measurement period and X% for the most recent 4 years. Most of the 
regression coefficients, including BI andBg are species specific. 
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This equation can readily be adapted for use in modifying an estimate of 
height growth to account for defoliation. Let PEHG be the proportion of 
expected height growth. If normal height growth is HGN, then height 
growth given defoliation is HGD = HGN(PEHG). Also, PEHG = HGD/HGN 
and In PEHG = In HGD - In HGN. 

Ferguson's equation can be used to compute normal growth and height 
growth given defoliation. That is, In HGN = I; BC it where X\ = X' 2 = 0, 
and In HGD = I ( - B& t where Xi > and^T 2 > 0. For a given tree, X' t = X t 
for i > 2. Therefore, 

In PEHG = (BiXi + 3*2 + . . . + B n X n ] - 



Imposing the simplifying assumption of letting .Xi -X^-X, we have 
PEHG 



(28) 



(29) 

The values for Si and B% are species specific. Their sums by species are 
as follows; for AF (also used for WF), -2.3661; for DF, -2.4757; for GF, 
-2.0008; and for ES, -2.9171. Figure 15 illustrates the curves. 

Equation 29 is used to compute the proportion of expected height growth 
on small trees. It is a 5-year model, which means that it is useful for modi- 
fying 5-year height growth. The defoliation measure is similar to that used 
in the topkill model except that the average is taken over 5-year intervals 
on the top third. 

Let PEHGhsT be the proportion of expected height growth on host h, tree 
size s(l = small in this case), at time T, when T is 5 and again when T is 
10. For a 5-year-long Prognosis Model cycle, 



exp (B h (1 - AVT hl5 }) 



(30) 
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Figure 15 Proportion of small-tree expected height 
growth plotted over defoliation proportion (from Ferguson 
1 988), The line for AF is also used for WF. 
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where B h are the species -sped fie constants described above, ati&AVTis 
computed in equation 23. If a 10-year-long cycle is used, the result of equa- 
tion 30 is divided by 2, the second 5-year part of the 10-year interval is 
computed, divided by 2, and the sum of these two estimates is the propor- 
tion of expected 10-year height growth. If a Prognosis Model cycle length is 
used that is not evenly divided into 5-year-long intervals, the model is ap- 
plied using as close to 5-year-long intervals as reasonably possible. 

PEHG haT is an average for the trees in a given class. The Budworm Mod- 
eling System assumes that the variances ofPEHG hsT are related to their 
averages like the variances of defoliation measures are related to their av- 
erages (equation 20). The available data, described in the section on using 
the Beta distribution, were used to verify this assumption. While the bases 
for making the assumption are not strong, the alternative of ignoring the 
variation between trees is a greater evil. Given the Prognosis Model esti- 
mate of height growth on tree n is HG' n , then 



H0 n = HG' n P>(PEHG hsT ,Q.Q3) (31) 

is the height growth on defoliated tree n, where tree n is of host h and size s. 

Medium and Large TreesK-Nichols (1988b) related growth to defolia- 
tion on large trees. His model for PEHG hs t is 

PEHGfot = 0.193 AV^f 814 -- 0212 ^> M l m ( 32 ) 

where Ml = PEHG hs(t _u when t > 1, and 1 otherwise. When s = 1 

AVT hlB (1,15 = i PRB hclt )/3T (33) 

for s = 2, AVT h2T = (I,!? = 4 pRB^iJ/ST, and for s 3, AVT h3T = ( = 7 
PRBhctiVST. Conditions controlling the summation on t described for equa- 
tion 23 also apply to equation 33. 

Equation 33 is computed once for each year of the Prognosis Model cycle. 
Figure 16 illustrates the function's behavior for several values of Ml. 
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Figure 16 Proportion of medium- and large-tree expected 
height growth plotted over defoliation proportion (from 
Nichols 1 988b), for various levels of M1 (see equation 32). 
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Diameter Growth 



Mortality 



The adjustment of the Prognosis Model growth estimate uses equation 31, 
that is, 

HG n = HG'n P (X* PEHG kst /T, 0.03} (34) 

Nichols' (1988b) model for proportion of expected diameter growth, which 
is measured in squared diameter increment, PEDDSh st > has the identical 
form as equation 32, However, the constants are different resulting in 
slightly different curves (fig. 17), The equation is 

= 0.084 A y^o.4725-o.0700Mi) M o.324i (35) 



where Ml = PEDDSh s (t-i) when t > 1 and 1 otherwise; and AVT hat is com- 
puted using equation 33. 

This equation is applied the same as the height growth equation (32). Di- 
ameter growth reductions are computed using this model for all three tree- 
size classes. 

The probability of individual tree mortality in budworm-infested forests 
has been studied with little general success. Both Ferguson (1988) and 
Nichols hoped to develop functions but received insufficient cooperation 
from nature. The budworms did not kill many sample trees, so the sample 
size was too small to support proposed analyses. Alfaro and others (1982) 
measured plenty of mortality on one stand near Pemberton, BC. They pub- 
lished two models that represent their data. The first relates mortality to 
tree size (PRMD, equation 36, fig. 18, line A), and the second relates mor- 
tality to cumulative defoliation (PRMC, equation 37, fig. 18, line B): 

PRMD n = 140.75 exp(-0.057.D n ) (36) 

PRMC n = 4.17 (10- n ) CMDfr 78 (37) 

where D n is d.b.h. in centimeters (for D n in inches, substitute -0. 1448 for 
-0.057), and CMD n is a 5-year cumulative percentage defoliation; both are 
for tree n, 
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Figure 17 Proportion of expected basal area Increment 
plotted over defoliation proportion (from Nichols 1 988b), 
for various levels of M1 (see equation 35). 
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Figure 18 Line A is probability of mortality plotted over 
d.b.h. (D, equation 36); and line B is probability of mortality 
plotted over cumulative percentage defoliation (CMD, 
equation 37). Lines A and B are from Alfaro and others 
(1 982). Line C is a refit of the data represented by line B; 
see equation 39. 



In a much more extensive study, Alfaro (1986) found that suppressed 
trees suffered significantly higher mortality rates than nonsuppressed 
trees. Slope and site quality were also related to mortality in budworm- 
defoliated stands. The number of years a stand suffered "severe" defolia- 
tion accounted for 47 percent of the variation in the mortality rate data. 
This finding is generally consistent with the model proposed from the 
Pemberton site, equation 37. That is, high values of CM) can only occur 
when there are several years of severe defoliation. 

In 1985, Crookston (1985, p. 216) published a combined version of equa- 
tions 36 and 37. This model has been criticized as being an unrealistic ex- 
trapolation of Alfaro and others' (1982) work. It predicts low mortality 
rates for large trees even when CMD is high. 

None of the models found in the literature are satisfactory for inclusion in 
the Budworm Modeling System. Equation 36, the function of diameter, as- 
sumes that stand structure and defoliation are the same as found at the 
Pemberton site. Equation 37 ignores site, competition, and tree-size effects 
that have been documented. Simple rates indexed by site or competition 
ignore the interactions with other variables. 

Clearly, an extensive study of budworm-caused mortality is needed. In 
the absence of such a study, and in the presence of work that has been ac- 
complished, the mortality model illustrated in figure 19 has been inserted 
into the Budworm Modeling System. The function used to represent the 
response surface is described below, after a biological interpretation is 
presented. 

RBAL n is the relative basal area in trees larger than tree n. It is com- 
puted by dividing the amount of basal area in trees larger than tree n 
(BAL n ) by the maximum basal area (BAMAX) that typically exists for the 
site: 
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Figure 19 Probability of mortality (PMRT] plotted over cumulative 
defoliation (CMD) and relative basal area in larger trees (RBAL} for 
several levels of d.b.h. (D). See equations 39-43. 



RBAL n ^BAL n /BAMAX 



(38) 



Values of RBAL n near zero indicate that tree n is growing with little com- 
petition. TVee n could be any size tree growing in a stand with little total 
basal area, or it could be the largest of many large trees. If RBAL n is near 
or above 1.0, then tree n is under competition. A small tree, growing under 
a large overstory in a stand that is at or above BAMAX, will have a high 
RBAL n , Site quality is inherent in this metric because a stand of moderate 
total basal area, growing on a poor site, could be close to its maximum. In 
the Prognosis Model, BAMAXia indexed by stand habitat type or it can be 
set by the model user. 

Figure 19 indicates that the probability of mortality, PMRT n , increases 
with increasing values of RBAL n . Increasing values of CMD n also increase 
the probability of mortality. For small trees, the mortality rate increases 
with increasing CMD n and increasing RBAL n faster than it does for large 
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trees. Furthermore, the maximum mortality rate for large trees growing 
without competition is lower than the maximum rate is for large trees 
growing with competition. 

The equation forPMRT n is an alteration of the Weibull p.d.f.; line C in 
figure 18 illustrates the shape of the function 

PMRT n = Bid - exp (- ( B2(CMD n + 1)) B3 (39) 

where El = 0.78, B2 = 0.003044, and S3 = 8.429. The values of B2 and B3 
were found by fitting Alfaro and others' (1982) published data points to 
equation 38, after arbitrarily setting the upper asymptote, Bl, to 0.78, 

Equation 39 is more flexible than equation 37. B2 controls at what point 
along the CMD axis the curve starts to climb, B3 is a power term that con- 
trols how steep the climb becomes before Bl flattens it off. 

In the final model Bl is a function of tree d.b.h. (D tl ). B2 and B3 were set 
equal to functions of D n and RBAL n that were crafted to provide the de- 
sired behavior. 

Let RBAL n = min (RBAL' n , 1), where RBAL' n is computed in equation 38, 
and let D n - min (D' n , 20), where D' n is the actual tree d.b.h. Then 

Bl n = 0.98- 0.01 D a ( 40 > 

XM n = 0.02 D n + (0.5 + 0.03 D n ) (I - RBAL n ) (41) 

B2 n = 0.00575 - 0.00645 XM n + 0.0031 XM$ (42) 

B3 n = 2.27- 2.94 XM n + 1,84 XM^ (43) 

Equation 39 is modified such that Bl, B2, and B3 are subscripted by tree n. 
The only remaining detail of the mortality model concerns how CMD n is 
computed. 

For each year in a Prognosis Model cycle, the overall average tree defolia- 
tion (CDEF hs T) is computed for host h, size class s, and year T, (Recall that 
Tis the number of years the model has run. For example, in the fourth 
year of a run, T = 4, in the tenth year, T = 10.) This is not defoliation of 
just the new foliage; it represents the defoliation in general, over the entire 
tree. It corresponds to the defoliation measure used by Alfaro and others 
(1982). ForT>5, ands = l 

CDEF hsT = 100 iL r-4 (1 - (I? = i la BAF hcat /? = : I a APB hcat }} (44) 

where APE and BAF are defined in equations 4 and 6, respectively. When 
s = 2, the summation over c includes 4, 5, and 6; and when s - 3, the sum- 
mation over c includes 7, 8, and 9. 
The maximum value of the CDEFhst is 

CMDF hs = max (CDEF hs5> CDEF hs6 ,..., CDJEF hsT ) (45) 

, CMDFh s has a theoretical range of to 500, where 500 means that a tree is 
completely stripped of all foliage each year of 5 consecutive years. The 
Beta distribution is used to generate an individual tree estimate; 

CMD n = 500 p (0.002 CMDF hs , 0.03) (46) 

where tree n is of host h and size s. 

Therefore, the mortality model used in the Budworm Modeling System 
includes elements of tree size, site quality, competition, stand size, and tree 
defoliation including variation in defoliation between individual trees in a 
tree-size class. The model estimate is assumed to be a 10-year mortality 
rate and is scaled (by discounting and compounding) for different Prognosis 
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Model cycle lengths. The Prognosis Model estimate of background mortal- 
ity is used in place of the Budworm Model estimate if the background rate 
is higher than the model rate. 

Damage Model The behavior of the damage model depends both on the damage model 

Behavior formulations and the behavior of the foliage dynamics model. By tracking 

the numerator and denominator for expressions of defoliation, the foliage 
dynamics model provides the most important variable controlling damage. 
The behavior of the foliage dynamics model has been illustrated for three 
defoliation sequences (cases 1-3, see table 11 and fig. 3-6). Recall that case 
1 portrays sustained heavy defoliation, case 2 portrays heavy defoliation 
that stops after 5 years, and case 3 portrays alternating light and heavy 
defoliation. All tree-size classes received the same initial defoliation profile 
in these illustrative runs, and in all cases the model was run for one 10- 
year-long cycle. 

Illustrations of damage model behavior for these three cases follow. The 
Bear Gulch stand data were run for each case, and outputs from the Bud- 
worm Modeling System's STATDATA option (Crookston and others 1990, 
p. 56) were used to generate the graphs. The only host species in this data 
set is Douglas-fir, so species differences described for some of the compo- 
nent models are not illustrated, 

Topkill The top part of figure 20 is a scatter plot of the probability of 
topkill plotted over tree height for individual trees. The symbol for each 
point corresponds to the defoliation case; that is, symbol 1 is for case 1 and 
so on. Note that small trees are those below 23 feet (7 m), and the prob- 
ability of topkill (equation 22) for them is different than for medium and 
large trees (equation 27). Therefore, it is not surprising to see different 
model behavior for trees below versus above the 23-foot threshold. 

Defoliation is measured as the maximum average over 5 years (equation 
24). The value of this measure is 0.95 for cases 1 and 2 and 0.65 for case 3. 
Given such a high average defoliation value and the shape of the Beta 
p.d.f. (fig. 10) when the mean approaches 0.0 or 1.0, the variation about a 
mean of 0.95 is quite low. That is, if the mean defoliation proportion is 
0.95, most of the individuals will fall near 0.95. The shape of the probabil- 
ity of topkill model for Douglas-fir (fig. 13) is very flat for the range of defo- 
liations 0.90-1.0. Therefore, the individual probabilities of topkill for cases 
1 and 2 are clustered on top of each other, forming a jumbled mass, near 
the upper asymptote of topkill function. For case 3, the probability of 
topkill is spread all over the plot, reflecting the variability in defoliation 
proportion when the mean is 0.65. 

For big trees, there is no tree-size effect on topkill probability (equation 
22 versus equation 27). On small trees there is an effect that can be no- 
ticed on trees below 23 feet in height. 

The bottom of figure 20 is a scatter plot of the number of feet killed over 
tree height for each defoliation case. As in the probability models, the 
amount of topkill models are different for small versus medium and large 
trees (equations 26 and 27). The linear effect for trees simulated under 
cases 1 and 2 implies that similar proportions of each tree's height are 
killed. The reason for this behavior is identical for the clustered behavior 
described for the probability of topkill model. For small trees, a distribu- 
tion function is used. Therefore, no correlation between proportion killed 
and any other variable is represented and should not be evident, 
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Figure 20 Probability of topkill (top) and amount of topklll (bottom) plotted 
over tree height for defoliation cases 1 through 3. There are nearly as many 
points for cases 1 and 2 as seen for case 3, but they plotted on top of each 
other forming a heaped mass. 
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Proportion of Expected Growth Trees that are not topkilled may 
sustain some degree of height growth reduction; all trees may sustain di- 
ameter growth reduction. Figure 21 illustrates the proportion of expected 
height and diameter growth reduction for each tree-size class and defolia- 
tion case. For diameter growth, the proportion is the same for all tree-size 
classes. This is an expected result because the same models are used, and 
the defoliation profile is identical, for all tree-size classes. Case 1 hits the 
trees the hardest, case 2 is second hardest, and case 3 is the lightest. 

The same general trend is true for proportion of height growth in me- 
dium and large trees. However, for small trees, case 2 results indicate that 
more of the expected height growth is retained than in case 3. This result 
is because (1) a different model is used for proportion of height growth in 
small trees (equation 29) than used for the others (equation 32), and (2) the 
small tree model takes the proportion in 5-year intervals. In case 2, defo- 
liation stopped after 5 years. Equation 32 models the recovery by includ- 
ing a lag variable that is missing in equation 29. Therefore, the odd behav- 
ior illustrated in figure 21 is an artifact of the model and not likely a 
reflection of nature. 

Scatter plots (fig. 22) of the proportion of expected growth plotted over 
tree size illustrate the effects of using the Beta distribution to represent 
variation between trees. As expected, there is overlap between cases, but 
trees run under case 1 are generally below those for case 2, and those for 
case 2 are below those for case 3. Note the reversal for the proportion of 
expected height growth on small trees; as expected, given the results dis- 
played in figure 21, the case 3 trees are below the case 2 trees. Also note 
that there are fewer points plotted on the height growth plot (fig. 22, bot- 
tom^ compared to the diameter growth plot (fig. 22, top). Trees that are 
topkilled are excluded from height growth reduction. 
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Figure 21 Proportion of expected diameter and 
height growth for each defoliation case by tree-size 
class. The odd behavior illustrated for small-tree 
height growth reductions Is due to an equation (30) 
being used on small tree? that differs from the 
equation (32) used for large trees. 
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Figure 22 Proportion of expected diameter (lop) and height (bottom) 
growth for each Prognosis Model sample free record by defoliation 
cases 1 through 3, There are many fewer trees plotted In the bottom 
graph because height growth reduction is relevant only on trees that 
were not topkllled. 
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Mortality Mortality is chiefly a function of cumulative defoliation. The 
5-year maximum CMD (equation 45) is charted for each tree-size class and 
defoliation case in figure 23. The most startling characteristics of figure 23 
are the high values of CMD for small trees compared to medium and large 
trees. A higher proportion of small trees' foliage is new compared to larger 
trees (table 9). Therefore, when constant proportions of new foliage are re- 
moved, the overall impact of defoliation is higher on small trees than on 
larger trees (review equations 44-46). The ranking between cases is as ex- 
pected: case 1 has the highest CMD, case 2 the second highest, and case 3 
the lowest. 

Figure 24 plots mortality rates over tree diameter for each Prognosis 
Model sample tree by defoliation case. The lower limits of mortality are 
those imposed by the Prognosis Model's own background rates. Case 1 
trees generally have the highest mortality rates. The rates are highly scat- 
tered due to the effects of the Beta, Small trees die at obviously much 
higher rates than larger trees. This is due to the high CMD values for 
small trees and the shape of the mortality function (fig. 19) that indicates 
that smaller trees die at higher rates than larger trees. 
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Figure 23 Cumulative percentage defoliation (CMD) 
for each defoliation case by tree-size class. Because 
smaller trees have a higher proportion of new foliage 
than large trees, values of CMD are higher for small 
tress even though percentage defoliation of new foliage 
is identical for all tree-size classes. 
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Figure 24 Probability of mortality plotted over tree diameter for 
each Prognosis ModGl sample tree record and by defoliation cases 
1 through 3. 



CONCLUSIONS AND RESEARCH NEEDS 

After reading the review draft of this manuscript, reviewer Michael 
Marsden (Rocky Mountain Station, USDA Forest Service, Fort Collins, CO) 
suggested being "more upbeat" in the conclusions, and offered the 
following: 

There ia a great deal of roam for improvement in the foliage and damage com- 
ponents of the Western Budworm Modeling System, The strong points of the 
present components are they are logical, simple yet complete, and can be easily 
modified. Some of the functions in the model may seem complex but they are great 
simplifications of the natural system. These functions are also selected to constrain 
the possible results to limits natural within the system. Adjustment can bo made 
to the model by adjusting coefficients as more information becomes available. 

Marsden's opinions about this model are as valid as mine. Much of the 
model is based on ideas and opinions. My original conclusions, after edit- 
ing, issue a stronger call for new work and care. They follow. 

The foliage and damage components of the Budworm Modeling System 
should be changed. The changes must result in behavior that better repre- 
sents what is known about tree and budworm interactions and dynamics; 
otherwise the model must stand. Research is needed to support creating 
the new ideas, analyses, and facts required to justify model modifications. 
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What can this model suggest about research priorities? The foliage dy- 
namics model is based on little data, yet I don't believe that research aimed 
at improving it is the highest priority. A sensitivity analysis that tests the 
overall model's performance relative to varying some of the relationships in 
the foliage dynamics model could provide convincing results to the con- 
trary. Clearly, the assumptions relating to defoliation measurements, and 
therefore all the ideas surrounding the adjusted potential biomass calcula- 
tions, directly impact the behavior of the damage models. 

Defoliation proportion is expressed as the amount of foliage on a tree af- 
ter defoliation, divided by the amount the observer expects there would 
have been if the tree had not been defoliated. Clearly, the divisor in this 
relationship is not directly observable. Yet the result, defoliation, is as- 
sumed to be field measurable without error in analyses relating defoliation 
to damage. Because most of the work on damage was done using some 
measure of defoliation, part of this model is necessarily devoted to predict- 
ing both the numerator and the denominator of defoliation relationships. 
These components may cause bias in the damage estimates, but I do not 
suspect that the qualitative behavior is wrong. 

In this model, height growth, topkill, and mortality are modeled sepa- 
rately. During the Canada/U,S. Budworms Program, some program par- 
ticipants (see the acknowledgments section) suggested that mortality is 
simply the ultimate form of topkilling. Furthermore, topkilling simply be- 
gins with 100 percent height growth loss. Any future analyses of height 
growth loss, topkilling, and mortality data might profit from this idea. In- 
cluding a model based on this idea in the present damage model is practi- 
cal. In any case, the topkill equations should be replaced with some based 
on analyses whenever they are ready. Building such relationships ranks 
high on a research priority list. 

A giant improvement over the modeling approach used here may require 
a totally different model. Photosynthesis creates the food for growth. 
Tracking photosynthetic ability with and without feeding would be a 
cleaner model form. What was done in this model involved attempts to 
track surrogate metrics of the fundamental processes. Using a physiologi- 
cal modeling approach would likely lead to unifying the foliage and damage 
components and lead to modeling topkill and tree mortality as a process. 

In the meantime, it is likely that this model will be used with few modifi- 
cations. A sensitivity analysis is needed both to demonstrate model behav- 
ior and component interactions and to provide guidelines to model users 
concerning how carefully they should measure some input data. Informa- 
tion regarding potential appropriate applications can be a result of sensi- 
tivity analysis. This research need ranks highest. 
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